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Abstract 

A theory is provided for the estimation of home ranges of animals from the standard 
mark-recapture technique in which data are collected by capturing, tagging and 
recapturing the animals. The theoretical tool used is the Fokker-Planck equation, its 
characteristic quantities being the diffusion constant which describes the motion of 
the animals, and the attractive potential which addresses their tendency to live near 
their burrows. The measurement technique is shown to correspond to the calculation 
of a certain kind of mean square displacement of the animals relevant to the specific 
probing window in space corresponding to the trapping region. The output of the 
theory is a sigmoid curve of the observable mean square displacement as a function 
of the ratio of distances characteristic of the home range and the trapping region, 
along with an explicit prescription to extract the home range form observations. 
Applications of the theory to rodent movement in Panama and New Mexico are 
pointed out. An analysis is given of the sensitivity of our theory to the choice of 
the confining potential via the use of various representative cases. A comparison 
is provided between home range size inferred from our method and from other 
procedures employed in the literature. Consequences of home range overlap are also 
discussed. 
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1 Introduction 



The use of space by animals is the result of a combination of internal factors, 
such as the physiology and morphology of the animal, and external factors 
such as the environment. It is well known that mammals, in order to conduct 
their daily activity, occupy on ly part of their available environment: the so- 
called home range (|BurtUl94.4 . A recent study on the scali ng of home range 



size as function of animal body mass or metabolic rate (|Jetz et all I2004TI 



shows that the home range dimensions are a compromise between two ingre- 
dients: the necessity for harvesting resources and the detection and response 
to intrusion. On one hand, the home range has to be large enough to meet 
energy requirements; on the other, it has to be small enough for the resident 
to be protected from intrusions of same-species foraging neighbours (Buskirk, 



2004). The importance of learning about home rang es stems not only from 



the intellectual need to understand animal movement (|Okubd . ll97i . but also 
from the practical value in the determination of the size of the home range: 
it is intimately related to a variety of ecological phenomena ran ging fr o m so- 
cial organization t o mating behavior and disease transmission (IWorfft 1 1997; 
Yates et all 120021 IParmenter and MacMahonl. Il983t Abramson and Kenkrd. 



20021: lAbramson et all l200i iKenkrel . 1200.4 120041 : iKenkre et all l2004 ) . There 



is already a body o f biological literature on home ranges and related animal 
movements (|Okubd . ll97i lMurravl . ll993| ). It is useful to build upon that liter- 
ature to provide a comprehensive mathematical description of the dynamics 
of ecological systems. 



Given the spatial probability distribution for an individual, the home range 
size has been typically defined as the contour that conta ins a fixed percentage 
(usually 95%) of the total volum e under the distribution ((Jennrich and Turner , 
1969; iFord and Krummel 1979i ). Home range size can be estimated from data 
of recorded locations over a sufficiently large period of time. These data can be 
obtained through a variety of methods including radio tracking and live trap- 
ping. Different techniques have been proposed for estimating home range si ze 
from location data of single animals (see e.g. the review bv IWortonl (|l987h . 
The various approaches can be divided into three categories. In the first, esti- 
mates are made using the peripheral points of the location data. In the second, 
the data are fitted to a pre-assumed probability distribution. In the third, the 
probability distribution is determined o nly from the sta tistical properties of 
the data (such as the one proposed by Anderson (1982)). The first method 
gives the maximum extent of the animal's range while the other two methods 
give a profile of the probability distribution inside the home range. 



In order to obtain reasonable accuracy, all three procedures have t he require- 

ment that the number of locations recorded for each animal be large ( Ford and Krummel 
19791 ) . Unless unusual efforts are made during the field sampling, the number 



2 



of loc at ions for each indi vidual is typically not very large (e.g. iMares et al. 
f)l980h : iBergstroml (|l988h ). On the other hand, the number of individuals 
is typically quite large. Home range estimati on from location data of many 
individuals, therefore, avoids such problems ( Ford and Krummel 1979f ): the 
positions are averaged over all the animals recorded. 



The purpose of the present paper is to develop a theoretical model which 
gives a simple prescription for the extraction of home range parameters from 
location data of animal population inside a limited region of space. Such a 
limited window of observation represents, for example, the size of the trapping 
array in a mark- recapture experiment, in which an animal is captured, tagged 



and then recorded every time it is recaptured ([Parmenter et all 200 



In the model we suggest, the motion of each animal is represented by diffusion 
in a confining potential, the latter representing the attraction of the animal 
to the home-place, the burrow. The potential has a characteristic width as- 
sociated with the size of the home range, which we call L. The underlying 
equation in our approach is the Fokker-Planck equat i on fo r the probability 
distribution for each individual (|OkuboL IT9791 iRiskenl . Il989i ). The stationary 
solution of this Fokker-Planck equation is used to calculate the infinite-time 
limit of the mean square displacement saturation value of all the individuals 
as function of L. Comparison with the measured mean square displacement 
allows then the determination of the home range size, expressed here in units 
of length. The home range size has been typically denoted in the literature 
in term of an area. The relation to our description is simply that that area 
is given by the product of L for the two directions. Application of our proce- 
du re to rodent measurements in Panama and in New Mexico may be found 



in ((Giuggioli et all l2004t lAbramson et all 120041 ) where our model in its sim- 



plest form has been successfully used to extract not only home range sizes but 
also diffusion constants from mark-recaptures of rodent populations. 



The practical output of our present theoretical procedure is a saturation curve 
for the observed mean square displacement as a function of L/G, the ratio of 
the home range to G, a length that is characteristic of the size of the observa- 
tion window. We predict a sigmoid shape for the saturation curve. An imme- 
diate consequence is that, for the greatest accuracy in the measurement of the 
home range, the observation window should be of the order of the home range. 
For certain potentials it is possible to write down simple analytical expressions 
for the saturation curve. For others the curve is obtained through numerical 
computation. Our theory also addresses the distribution of home ranges ac- 
cording to habitat, equivalently home range overl ap, a quantity ind ependently 



accessible through allometric scaling arguments ([Jetz et all 12004). 



The paper is organized as follows. The general problem of calculating the aver- 
age mean square displacement for a population of individuals, each one living 
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in its own home range, and observed only inside a spatially limited window, is 
addressed in Sec. 2. The sensitivity of the saturation curve to the choice of the 
confining potential is studied in Sec. 3 through various representative cases. 
The case of a non-uniform distribution of home ranges and considerations for 
experimentally determining the average inter-home distance (related to the 
home range overlap) is the subject of Sec. 4. The comparison between home 
range size inferred from our method and the so-called convex polygon method, 
usually employed in the literature, forms Sec. 5, and conclusions are in Sec. 6. 



2 Mean square displacement in a probing window: general consid- 
erations 



We model the motion of an animal living in its home range through the Fokker- 
Planck equation for the probability distribution V(x, t) 



dV(x, t) 
dt 



d_ 

dx 



dU{x) 
dx 



V(x,t) 



+ D 



d 2 V{x,t) 
dx 2 



wherein D is the diffusion coefficient of the animal and U (x) is the potential 
in which the animal is forced to roam. The potential U (x) is a representation 
of the bias or reduced randomness associated with the walk. A pure random 
walk as in a simple diffusive process has U(x) = 0. When U(x) ^ we identify 
its characteristic length with the home range size L. 



The 2-dimensional counterpart of Eq. (1), written in polar coordinates, is 



dV(r,(f),t) 
dt 
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(2) 



The description provided by (2) would be appropriate for a wide range of 
animal motion contexts. The third dimension is very rarely required but can be 
easily incorporated. All the essential concepts are, however, easily represented 
through the description provided by the 1-dimensional version. Therefore, we 
restrict ourselves in the present paper to Eq. (1). It is straightforward to 
generalize all considerations to higher dimensions if required. 

The overall characteristics of the motion can be obtained by calculating just 
the moments of the distribution V(x,t) rather than the full V(x,t) which 
entails an integration of the probability. In a typical experiment a probing 
window is used, i.e., the animal is observed only inside a limited region of 
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space, as in a trapping array in a mark-recapture experiment. The integration 
to calculate the moments is then performed only over the probing window. 
The second moment of V(x, t), i.e., the mean square displacement, is given by 



Ax 2 {t) 



I—G/2 dx {x 



xo) 2 V Xo (x,t) 



I—G/2 dx "Pxq {x, t) 



(3) 



where G is the dimension of the window and Xq is the position of the animal 
at time t = 0. Because initially each animal can be anywhere inside G, the 
numerator and denominator of Eq. (3) have to be averaged over all the possible 
initial positions inside the window. We then have, for the average, 



Ax 2 {t) 



I-G/2 dx I-G/2 dx ( X ~ X o) 2 Px^xoix, 0) V Xc 

I—G/2 dx ® I—G/2 dx 'PxcXq {x, 0) "PxcXq {x, t) 



(x,t) 



(4) 



We have introduced here the label x c to represent the burrow position of each 
animal. Equation (4) is the contribution to the mean square displacement 
of an animal whose burrow is at x c . A further average over the distribution 
of burrow positions is necessary. If this distribution is denoted by p(x c ), the 
observed mean square displacement within the window of size G is given by 



Ax 2 (t) 



I-oo dx c p(x c ) I_Qj2 dXQ I_g/2 dx ( X 



•^o) ^Px c ,xo{Xi^)l^x c ,xo^X^t) 



I-oo dx c p(x c ) I_Gj2 dx$ I-G/2 dx T^XcXo { x i 0) ^Px c ,x {x, t) 



r G/2 



(5) 



Short -time measurements of Ax 2 (t) can be used (jGiuggioli et al!l2004tlAbramson et al 
2004) to obtain the diffusion constant D. In the present paper we are inter- 
ested only in the home ranges, consequently in the infinite time limit of Eq. 
(5) which requires only the steady state solution. 

Analytic solutions of the Fokker-Planck equation for all times are known only 
for very few cases of U(x). However, steady s tate sol u tions are known for any 
potential explicitly in t erms of an integral (Risken, 19891 : Kus and Kenkrd . 
1992HParris et al! l200lh 



V Xc , X0 (x,t -> +00) 



-U(x-x c )/D 



Equation (5) for t 



JZodx'e-^')/ ' 
00 can, thus, be written as 

r G/2 7 rG/2 , / \ o U(x -x c ) + U(x-x c ) 



(6) 



I-oo dx c p(x c ) I_ G/2 dx I_ G/2 dx (x - x y 



I-oo dx c p(x c ) I_q/2 d x I-G/2 dx e 



U(xq—x c )-^U{x — x c ) 
D 



(7) 



and further reexpressed in terms of quantities related to moments of the steady 
state probability density, equivalently of exp [— U(x)/D]: 
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Axl = 2< 
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The expression (7) can be reduced further if the burrow distribution p(x c ) is 
uniform in space. We obtain 



G J G G dy y 2 g(y) - (j G dy y 3 g(y) - J° G dy y*g(y) 
G I G G dyg(y) - (j G dy yg(y) - J° G dy yg(y) 

where g(y) is the convolution of exp [— U(x)/D] with itself: 

g(y) = 



(9) 



U[x c ] + U[x c -y] 

dx r e D 



(10) 



If the potential U(x) remains finite for all finite values of x, g(y) = g(—y). 
and it is possible to write a simpler form of (9): 



Jo G dy(G - y)y 2 g(y) 



J G dy(G-y)g(y) 



Clearly, the trapping array of width G can be interpreted as a probe into 
the system whose characteristic width is L. The mean square displacement 
depends on the relative magnitude of L and G of the probe. In the limit of an 

infinitely large probe, ( Ax 2 J measures simply the characteristic length of 
the system. It is thus natural to define the home range length L for arbitrary 
potentials as the square root of the limit G — > oo of Eq. (9), 



J^dyy*g(y) 
\ f^dyg(y) ' 



(12) 



when such a limit exists. If the probe is very small compared to the home 
range width, the trapping array will measure a quantity associated with the 
width of the grid. In fact g(y) becomes a constant in the limit L — > oo and 
Eq. (11) gives G 2 /6 for the mean square displacement . 



In a recent article, one of the present authors ( Kenkrel . 20051 ) has given an 
alternate formulation in terms of the Fourier transform of the steady state 
probability distribution (6). It has been shown there that the mean square 
displacement is given simply in terms of Fourier-space integrals of the product 



* While the definition (12) of the ho me range L is the mos t natural, alternate 
definitions are possible as used, e.g., in ( Abramson et al. . 20041 ). 
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of the square of the sine function (which carries information about the probe) 
with respectively the square, and the derivative of the square, of the transform 
of the steady state probabilities (which carry information about the home 
ranges). Alte rnative expres sions equivalent to Eq. (9) and (12), respectively, 
are given in ( Kenkrel . 2005) as 



Ax 2 



S-oo dk 



d 2 P 2 (k) [l-cos(Gfc)] 
dk 2 A? 
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L 
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O io dyy z e d 

V(y) 



k=o /o dye n 
where P(k) is the Fourier transform of exp [—U(x)/D]. 



(13) 



(14) 



Clearly, a dimensionless quantity of crucial importance to the analysis is the 
ratio ( of the home range to the observational probe length G: 



C = L/G. 



(15) 



In the next section we study the functional dependence of the saturation curve 
on this quantity ( with attention to the effects of the details of the confining 
potential, assuming that p(x c ) is a constant. 



3 Dependence on the details of the confining potential 



The precise shape of the confining potential U(x) obviously depends on the 
detail of animal movement, such as habitat and distance between neighbours. 
Since this detail is largely unavailable, it is important to determine the sensi- 
tivity of the deduced value of the home range size L to the choice of U(x). It is 
clear that, when plotted as a function of ( (see Eq. (15)), the mean square dis- 
placement Ax 2 s for each potential starts out as L 2 when £ — > and saturates 
to G 2 /6 when ( — ► oo. 

Extensive studies we have carried out with different potentials have made 
it clear that the curvature at the bottom and the steepness with which U(x) 
becomes infinite both play a role in shaping the saturation curve of Ax 2 s . More 
precisely, the rise of Ax 2 s is controlled by the steepness of the potential when 
i/I>1: the steeper the rise to infinity of the potential, the smaller the value 
of ( for which the saturation curve grows faster than the L 2 dependence at 
L = 0. In addition, the curvature of U(x) for i/L <C 1 determines the way Ax 2 s 
approaches the value G 2 /6: the larger the curvature, the slower it approaches 
the asymptote G 2 /6. In Fig. 1 we show this dependence by comparing four 
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characteristic potentials: a box potential, a harmonic potential and two types 
of logarithmic potentials. 




Fig. 1. Mean square displacement at saturation for four different potentials: a box 
potential (dashed line), a harmonic potential (solid line), and two logarithmic po- 
tentials (see text for definition) one with n = 3 (dotted line) and the other with 
n = 8 (dash-dotted line). The inset shows the corresponding stationary distributions 
V{x) as obtained from Eq. (6) by putting x c = 0. 



3. 1 Box and harmonic potentials 



The box potential has the steepest rise among the four curves in Fig. 1, since it 
diverges at a finite distance. Observe also that it is the slowest to reach G 2 /6, 
being the one with the smallest curvature (zero) clo se to the origin . Such 



a potential has been c onsidered in our previous work ( Giuggioli et al. , 120041: 



Abramson et all 120041 ) for extracting home range sizes from mark-recapture 



data for two different rodent populations. In those investigations the satu- 
ration curve was numerically simulated. However, it is possible to calculate 
analytically all the integrals in Eq. (8) by selecting judiciously the limit of 
integration for the variable x c as function of the relative dimension of the box 
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width and the probe width. The resulting expression in terms of ( is given by 



G 2 /Q 



18C 2 (5-3^) 

5 (3-Vec) 

3 (3-5%/6Q 
5 (l-3^<) ' 



for C < VQ, 
for C > V6- 



(16) 



The limiting behaviour is that Axl s 
4/(15^0 for large (. 



L for small ( and Ax 2 



G 2 /6 



The h armonic potential has been used in one of our previous studies ([Abramson et al 
2004) for extracting the home range size for the deer mouse, Peromyscus 
maniculatus, in New Mexico. The shape of the potential is given by U(x) = 
D(x/L) 2 . Since the domain where U(x) is not zero extends over the en- 
tire real axis, it is convenient to calculate the convolution (10) which gives 
g(y) = L^Jtc/2 exp [— y 2 / (2L 2 )}. Integrating Eq. (11) gives the mean square 
displacement at saturation as 



G 2 /Q 



6C 



sinh 


a/20 2 ; 




sinh 


(1/2C) 2 


-iv^ e(1/2C)2erf (v%) 



(17) 



with Ax 2 s ~ L 2 when ( < 1 and Ax 2 s ~ G 2 /6 [1 - (7/120C 2 )] when ( > 1 . 
The intersection in Fig. 1 between the harmonic potential curve (solid) and 
the box potential curve (dashed) becomes evident here given that the former 
approaches G 2 /6 linearly while the latter approaches it quadratically as £ — > 
oo. 



3.2 Logarithmic potentials 



As an example of U(x), whose corresponding steady state distribution given 
by (6) decays to zero slower (algebraically) than in the box or the harmonic 
case, we consider a family of potentials of the form 

{2 ~\ n 
1 + -^7T2 , (18) 
(«n£) J 

with K n = — 3/2 and n > 2. These potentials have a quadratic dependence 
when i/L < 1 and a logarithmic dependence for x/L ^> 1 such that V(x) 
~ x~ 2n as |x| — > oo. It is straightforward to check that the expression (18) 
is consistent with the definition (12) of the home range L for any value of n. 
In the limit of n — >• +oo (18) reduces to the harmonic case. We have studied, 
in particular, the cases from n = 2 to n = 8 and have obtained analytical 
expressions for the mean square displacement. We do not display them here 
because they do not add to the understanding. In Fig. 1 we show the cases 
n = 3 and n = 8. Already for n = 8 the saturation curve for the logarithmic 
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0.0 0.5 1.0 1.5 

o/G 



Fig. 2. Mean square displacement at saturation for the potential 
U(x) = Dln(l + (x/a) 2 ). Compared to the potentials depicted in Fig. 1, 
here U(x) grows to infinity qualitatively slower. The corresponding distribution 
V(x) does not possesses a finite second moment. The long-tailed distribution 
changes drastically the behaviour of Ax 2 s when very large grid sizes are used. The 
growth of the saturation curve for small a/G is linear and not quadratic. 

and the harmonic cases are very close to each other. The inset of Fig. 1 shows 
the origin of such similarities by comparing the corresponding steady state 
probability distributions V(x). 

As mentioned above, the V(x) associated with the logarithmic potentials pos- 
sesses long tails. Probability distributions with long tails can be appropriate 
when the motion of the animal cannot be represented by a simple random 
walk. More complex types of walk may occur if the walker awaits for very long 
times between jumps, or if the jumps are of very large distance. Long-tailed 
V(x) are characterized by the feature that certain moments of the distribution 
become infinite. If all the moments beyond the first are infinite, Eq. (12) is no 
longer applicable for defining the home range width. A qualitatively different 
behaviour is expected for large values of G. To illustrate this situation, we 
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consider the potential 



,2 

U{x) =L>ln(l + ^-) (19) 



whose corresponding V(x) is the Cauchy distribution ( Lemonsl . l2002j ) V 



[1 + (x/a) 2 ) Also for this case the mean square displacement can be 
obtained analytically and it is given by 



At 2 
Xss 6 



G 2 /6 



^ -4£ 2 



tan-Mi) _^i n (i + 



(20) 



where f = er/G. As ^ -> 0, Ax 2 J(G 2 /6) ~ (12/ir)a/G in Eq. (20). A linear 
growth of the saturation curve emerges, as depicted in Fig. 2. Notice that L, 
as defined by Eq. (12), does not exist for this potential. The home range is 
therefore defined as the characteristic length a, the ratio £ being the coun- 
terpart of ( of (15). This different qualitative behaviour with respect to the 
previously analyzed cases could be exploited for determining the character- 
istics of the animal walks by making various measurements with large grid 
size. A sufficient number of these measurements would allow one to discern if 
the saturation curve is growing quadratically (as in the examples of Fig. 1) or 
linearly (as in Fig. 2). 

The examples of this section illustrate that the choice of U(x) in any applica- 
tion of the present theory should be assessed in each case, based, for example, 
on a priori knowledge of the specific animal behavior. The different potentials 
shown here give an overview of the possible qualitative behaviours of the mean 
square displacement at saturation. 



4 Inhomogeneous distribution of home ranges: the case of a peri- 
odic arrangement 



The results obtained in Sec. 2, from Eq. (9) onward, and the examples devel- 
oped in Sec. 3, assume a continuous and homogeneous distribution of burrow 
location x c . A more realistic situation invokes the home ranges arranged in a 
non-continuous manner, the centers of adjacent ranges (the burrow locations) 
separated by some characteristic distance a. In this section we show an exam- 
ple of how this feature may be incorporated in our analysis and and how, in 
principle, a can be deduced from mark-recapture measurements. 

Let us suppose for simplicity that the home ranges are distributed in a peri- 
odic array, with a being the distance between nearest neighbors. The mean 
square displacement, measured within a window of linear size G, is now a 
function of a and L, besides depending parametrically on G. As in Section 
2, G can be used to rescale the two coordinates a and L. The function 
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L/G 



Fig. 3. Contour plot of the normalized mean square displacement, Ax% s /(G 2 /6), as a 
function of the normalized home range size, L/G and the normalized inter-home dis- 
tance, a/G. The lines are the result of simulations of the harmonic model (Gaussian 
occupation of space). The fluctuations in the lines are an artifact of the construction 
of the contours from discrete simulations. 

Ax 2 s /(G 2 /Q) = f (L/G, a/G) is universal and does not depend on the size 
of the observation window. We show a contour plot of this function in Fig. 3, 
as calculated by numerical simulation of the harmonic model (Gaussian prob- 
ability distributions). At a = 0, the shape of the surface coincides with the 
curve calculated in Eq. (17). The contours of equal mean square displace- 
ment are nearly vertical lines in this plot, indicating that the dependence on 
the inter-home distance a is very weak (in particular for small values of the 
normalized mean square displacement), except for a region well defined in a, 
where the contours shift from one value of L to another. This feature will cer- 
tainly be of relevance if an experiment is designed in order to measure both a 
and L, since the uncertainty on a will tend to be large. 

In general, given that the function Ax 2 s is nonlinear in both its variables a/G 
and L/G, two or more measurements are necessary to determine the home 
range size and the inter-home distance. In Fig. 4 we show a hypothetical 
situation in which three measurements are supposed to be taken on the same 
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0.8 



1.0 



L 



Fig. 4. Graphical explanation of the procedure to determine home range overlap 
via a hypothetical measurement. Both L and a can be determined by taking three 
measurements of the same population, using three windows, of sizes G±, G 2 and 
G3. The measured mean square displacements are supposed to be Azj = 0.045, 
Ax| = 0.029, and Ax^ = 0.041. The corresponding contours intersect at L = 0.5, 
a = 0.375, providing these values as a result. Both a and L are displayed in units 
of G (arbitrary linear units). 

population, using three windows sizes, G± — 1, G 2 — G\/2 and G3 = 3Gi/4. 
The results of the measurement are curves of constant Ax^ s in the plane (L, a). 
With three measurements, Arrf , Ax?; and Ax%, three curves are obtained, and 
the model predicts: 



This is a system of three (nonlinear) equations with two unknowns, L and a, 
and its solution can be found as the intersection of three curves. These curves 
are shown in Fig. 4, displaying an intersection very near L = 0.5, a = 0.375 
(within the accuracy of the fluctuations of the contours). The curves were 
obtained from the normalized function, shown in Fig. 3, using the appropri- 
ate contours. As mentioned above, the weak dependence on a may hinder its 



Ax\/{Gl/Q) 



f(L/G 1 ,a/G 1 ), 
f(L/G 2 ,a/G 2 ), 



Axl/(G 2 2 /6) 




(21) 
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a = a/L 



overlap = 50 % overlap = % 

B ^ C g ^ C 

" A " " A * 

a = 3L/4 a-L 
a = 3/4 a = 1 



Fig. 5. Illustration of home range overlap situations, and the variables characterizing 
them. The schematic diagram shows three cases: 100%, 50%, and 0% overlap, with 
arrows representing the extent of the home range of three neighbouring animals, A, 
B and C. The plot shows the relation between the overlap, the exclusivity of space 
use, and the normalized inter- home distance a = a/L. 

determination by the present method. It is clear that the choice of the ap- 
propriate values of the window sizes is critical to obtain the best results. This 
must be done specifically for each situation, with the help of an informed guess 
of the range where both a and L lie. Regardless of this practical difficulty, the 
procedure we describe provides a method for an immediate measurement of 
an important quantity that is hard to obtain by other means. Additionally, 
if the population is not well characterized by a typical inter-home distance, 
as we suppose here (for example, if the inter-home distance is bimodal due 
to gender differences or other polymorphism), the model can be immediately 
modified to incorporate those features. 



Moreover, the inter-home distance is closely related to the overlap of home 
ranges (or to the exclusivity of space use). See Fig. 5 for an illustration of 
three typical situations. The home ranges of three neighbouring animals, A, 
B, and C, are displayed as arrows, schematically representing t he extent of th e 
area occupied by 95% of the norm of V(x), as usually defined ([WortonL 11987). 
When a = L/2, the first neighbours, B and C, of A have their home ranges 
situated at the border of A's home range. Then, A does not have exclusive 
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use of any part of its own home range (overlap equal 100%). On the other 
extreme, when a = L, the home ranges of B and C are completely outside 
A's. In consequence, the exclusivity of A is 100% (home range overlap 0%). 
This value of exclusivity, certainly, is maintained for any value of a > L. An 
intermediate situation, in which the exclusivity of animal A is equal to 50% 
(as well as its overlap with the neighbours) is also shown. 



The exclusivity of space use has recently b een found t o obe y an allomet- 
ric scaling relation with the animal mass bv iJetz et al.l ([20041 ) . The present 
theory provides a method to determine both the home range size and the 
home ra nge overlap, and thus to verify the scal ing of these and related mag- 



nitudes ( West et al. . 1997 : Banavar et al. . 1999). 



5 Comparison with convex polygon calculations 



Home range sizes have been often deduced from the measurement of the so- 
called minimum convex polygon of an ani mal position. A lthough this proce- 



dure suffers from a number of drawbacks ([Wortonl . 119871 ). it is used rather 
widely. Among its drawbacks it is easy to recognize at least a logical and a 
methodological one. The surveyed perimeter provides no information about 
the use of space inside it, effectively encompassing areas that may be in- 
accessible to the animal, or potentially huge areas of very low frequency of 
utilization. The methodological one is the fact that the measured area con- 
verges very slowly to the actual home range, and as such the observation of a 
few tens of positions provides a very bad esti mation. B o th kin ds of flaws have 



been recognized in the literature before (see IWortonl ([19871) and references 
therein), and suggestions have been made to compensate for them. These pro- 
posals, such as discarding some fraction of the extreme positions from a set 
of observations to compensate for the first, or to join the perimetral points 
in a fashion different from the minimum convex polygon, are surely arbitrary. 
Furthermore, methodologically, they are obviously subject to uncontrollable 
errors. 

In the following we illustrate how the calculation of the mean square dis- 
placement we have given in the present paper provides a rapidly converging 
measurement of the home range size. This is an additional advantage when 
displacements of animals belonging to a population are more readily accessi- 
ble than repeated measurements of the position of individual animals. For the 
sake of the illustration, we will consider that the probability of space use of 
an animal is a bivariate Gaussian distribution of variance a: 

1 * 2 +y 2 

V ^ y) = 2^ 6 ^- (22) 
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We consider this symmetric distribution for simplicity, but the discussion ap- 
plies equally to a distribution with anisotropic a. Moreover, the conclusions 
are equivalent for more general distributions, including those with a finite 
cut-off. 

If we define the home range of the animal whose space use distribution i s 



described by (22) as the area A that contains the 95% of V (jWortonl [1987), 
a simple integration gives 

A = ttR 2 = 2trx 2 (- In 0.05) « 18.8a 2 . (23) 

This the quantity that we intend to measure by both methods. In Fig. 6 we 
show the area of the minimum convex polygon defined by a set of N points 



t 1 i 1 i 1 r 




Fig. 6. Area A enclosed by the minimum convex polygon corresponding to a set of 
random points with Gaussian distribution in the plane, as a function of iV,the num- 
ber of points of the set. The line shows the average of 50 independent realizations. 
The inset shows a detail of the same function, with arrows showing the necessary 
number of points to obtain a good measurement of the home range, defined as the 
95% of the space occupation which, for a = 1, and A 18.8, is N = 135 points. 
The variance a has arbitrary units of length, and A those of length squared. 
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drawn at random with a bivariate Gaussian distribution of o = 1. Observe, 
firstly, that the area grows unboundedly as the number of observations grows, 
since V is unbounded. More relevant from the practical point of view is the 
fact that the growth is very slowly, and that the area A (obtained with 95% 
rule and marked in Fig. 6 with an arrow) is achieved after the observation 
contains, on average, 135 points (see the inset of Fig. 6). 

Now consider that one wants to determine A by measuring displacements, 
instead of positions. The relevant quantity to be estimated is the variance 
a, immediately derived from the home range size L, accessible through our 
present theory explained in Sec. 2. Expression (23) gives then the home range 
area. Figure 7 shows that the measurement of the variance, <jn, for a finite 
set of N observations, converges very rapidly to the actual value (which is 1 
in this case), when the number of observations is increased. Indeed, with just 
10 observations of the displacement the variance can be estimated with an 
accuracy greater than 95%. Ten displacements correspond to just 11 positions 
of a single animal, if taken at intervals long enough that they are uncorrelated, 



1.00 




N 



Fig. 7. Measurement of the variance of a bivariate Gaussian distribution (a = 1), as 
a function of the number of displacements used for the estimate. The curve shows 
an average based on 104 realizations. 
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or 10 displacements of different animals for which an average distribution 
would be found. 

In summary, the mean square displacement provides a faster convergence to 
the area of the home range than the construction of the convex polygon. In 
addition, the estimation of the distribution V (unimodal in this case, but easily 
generalizable) by that method provides information about the use of the home 
range, which is inaccessible to the convex polygon procedure. 



6 Conclusions 

The determination of home range dimensions and spatial overlap of two neigh- 
bouring home ranges from field observation is a subject of great interest for 
the understanding of animal motion. We have provided here a general theory 
to extract such demographic parameters on the basis of the measurement of 
displacements of individual animals in a population. 

The most common techniques for gathering information of home range size 
employ trapping of animal and radiotelemetry observations. Our theory in 
the present paper has been constructed with the specific goal of interpret- 
ing data obtained from the former type of observation, i.e., mark-recapture 
experiments, in which the sampling area is finite. 

The motion of an animal inside its own home range has been modeled by 
a Fokker-Planck equation (1), i.e. by diffusion in a confining potential U(x). 
While the equations considered are, for simplicity, 1-dimensional, extension to 
higher dimensions is straightforward and unnecessary for practical purposes. 

Even though the precise determination of the home range size L depends 
on the choice of the potential U(x), the general sigmoid shape of the satu- 
ration curve in our theory indicates that the difference in the results is not 
substantial if the window size is chosen such that L < G. Eventually, the 
choice of the right potential is to be determined for each given case, on the 
basis of biological information of the animal population under study. We have 
shown that, for those situations in which the second moment of the distri- 
bution V(x) = exp(—U(x)/D) is not finite, the saturation curve for large G 
grows linearly with L rather than quadratically as in more conventional po- 
tentials (box or harmonic). Such cases may arise when the animal motion is 
not simple but involve a more complicated random walk such as a Levi walk 
or flight. This means that our theory of mark-recapture observation may be 
used, in principle, to determine whether the animal population is performing 
a Gaussian random walk or a more complicated walk. By measuring the mean 
square displacement at saturation with different values of the probe length G 
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(sufficiently lager than L) it might be possible to determine if the saturation 
curve grows quadratically or linearly. 

The obtained parameters that characterize the average use of space, when 
obtained via mark-recapture observations and their interpretation with our 
present theory, converge rapidly to the expected values. We have shown in 
Sec. 5 that this is not the case in the application of the traditional minimum 
convex polygon method. We thus suggest that displacement measurements 
of an animal population should be considered as the most appropriate way 
for determining home range dimensions if radiotelemetry methods are not 
available. 

The other important demographic parameter that can be extracted from mark- 
recapture measurements is the inter-home characteristic distance of the an- 
imals, called a in the present paper. Such a length is simply related to the 
mean overlap (or the mean exclusivity of space usage). We have outlined a 
procedure to extract this parameter quantitatively from mark-recapture ob- 
servations and provided a general way to verify directly from mark-recapture 
experiments the scaling of home range overlap as function of body mass. 
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